Interplay between parallel and diagonal electronic nematic phases in interacting 

systems 



O 
O 

(N 

S3 

a 



Hyeonjin Doh,* Nir Friedman, and Hae- Young Kee^ 
Department of Physics, University of Toronto, Toronto, Ontario M5S 1A7, Canada 

(Dated: February 6, 2008) 

An electronic nematic phase can be classified by a spontaneously broken discrete rotational sym- 
metry of a host lattice. In a square lattice, there are two distinct nematic phases. The parallel 
nematic phase breaks x and y symmetry, while the diagonal nematic phase breaks the diagonal 
(x + y) and anti-diagonal (x — y) symmetry. We investigate the interplay between the parallel and 
diagonal nematic orders using mean field theory. We found that the nematic phases compete with 
each other, while they coexist in a finite window of parameter space. The quantum critical point 
between the diagonal nematic and isotropic phases exists, and its location in a phase diagram de- 
pends on the topology of the Fermi surface. We discuss the implication of our results in the context 
of neutron scattering and Raman spectroscopy measurements on I^-ajSr^CuCU. 

PACS numbers: 71.10.Hf,71.27.+a 
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I. INTRODUCTION 

Recently, there has been a great effort to understand 
intrinsic phases of a doped Mott insulator in the con- 
text of high temperature superconductors. It has been 
proposed that quantum fluctuations of a Mott insula- 
tor introduced by hole doping lead to intermediate forms 
of matter, dubbed as electronic smectic and nematic 
phases. 1,2 In analogy to classical liquid crystals, the smec- 
tic phase breaks translational symmetry along one direc- 
tion, while the nematic phase breaks rotational symme- 
try. 

The evidence of such inhomogeneous and/or 
anisotropic liquid phases has been found in strongly 
correlated electron systems. 3-8 In particular, the clear 
evidence of a nematic liquid phase has been reported 
in two-dimensional electron gases in magnetic fields in 
ultra-clean samples. 5 The observed strong anisotropy of 
longitudinal resistivity has been explained by the onset 
of a nematic phase at low temperatures. A recent theo- 
retical study of the nematic phase using a quadrupolar 
interaction, F%, has offered non- Fermi liquid behavior in 
the nematic phase as well as near the quantum critical 
point. 9 This is originated from large fluctuations of the 
overdamped collective modes within the RPA theory. 
A non-pcrturbative approach using higher dimensional 
bosonization reproduced the quantum critical behavior 
with the dynamical exponent of z = 3, and verified the 
non- Fermi liquid behavior in the nematic phase. 10 

However, it was shown that the nature of phase tran- 
sition of the model with the quadrupolar interaction is 
quite different when we take into account an underlying 
square lattice. 11,12 On a lattice, a nematic phase can be 
achieved via a spontaneously broken point-group symme- 
try due to interactions between electrons. For example, 
it can break x and y symmetry of a square lattice. An 
essential consequence of nematic order is a deformation 
of a Fermi surface. It was shown that the transition from 
isotropic liquid to the nematic phase which breaks x and 
y symmetry of the square lattice is strongly first order at 



low temperatures. The nematic order parameter jumps 
at the transition to avoid the van Hove singularity, thus 
suppressing the Lifshitz transition. The transition takes 
place at arbitrarily small attractive quadrupolar interac- 
tion at the van Hove band filling. The transition changes 
to a continuous one at a finite temperature, but is not af- 
fected by either the next neighbor hopping, t', nor small 
dispersion in the third direction. 

In this paper, we study two distinct nematic phases in 
the square lattice, and investigate the interplay between 
them. The parallel nematic phase previously studied 12 
breaks a symmetry between x and y, while the diagonal 
nematic phase breaks a symmetry between two diagonal, 
(x + y) and (x — y) directions. The order parameter as- 
sociated with the parallel (A) and diagonal (A') nematic 
phases are defined as follows: 



A = F 2 y^(cos k x - cos k y ) ( 4.c k ) , 

k 

A' = i^^^sin^sinfcy) (c k c k \ . (1) 



It was shown that A' is always for a given quadrupolar 
interaction F2, which implies that a preferred direction 
for electron momenta has been selected to be parallel 
to the crystal axes. Here we study the interplay between 
parallel and diagonal nematic phases using a phenomeno- 
logical model with two different strengths of interactions, 
F2 and F' 2 for the parallel and diagonal nematic orders, 
respectively. We find that the transition to the diagonal 
nematic ordered state from isotropic liquid phase occurs 
above a critical value of interaction F%, and it is second 
order as a function of chemical potential. The compe- 
tition between the diagonal and parallel nematic phases 
leads to suppression of the strengths of both phases, while 
they coexist in a finite window of chemical potential. 

The paper is organized as follows. We describe the 
effective model Hamiltonian for the nematic order in sec- 
tion II. The mean field analysis at zero temperature is 
given in section III. The effect of t' is also presented. 
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Wc discuss the implication of our results and compare 
with neutron scattering and Raman spectroscopy mea- 
surements on La 2 -a;Sr 2 ;Cu04 in section IV. We provide 
the summary of our findings and future works in the last 
section. 

II. HAMILTONIAN FOR NEMATIC ORDERS 

Within a weak-coupling theory, the instability toward 
a Fermi surface deformation, often referred to as Pomcr- 
anchuk instability, 13 has been discussed in a Fermi liquid, 
t — J model, Hubbard model, and the extended Hubbard 
model. 14-18 It was shown that a strongly nematic phase 
(quasi-lD) is stable in a strong coupling limit of the two 
dimensional Emery model. 19 A phcnomenological model 
with quadrupolar density interaction in the continuum 
case was introduced in Rcf.9, and it was extended to the 
square lattice in Ref.11,20. 

The quadrupolar density interaction involves two dis- 
tinct nematic phases in the square lattice. For the par- 
allel nematic phase, the Fermi surface expands along 
k x (k y )-axis and shrinks along the k y (k x )-axia. On the 
other hand, for the diagonal nematic phase, the Fermi 
surface expands along (k x + k y ) and shrinks along the 
(k x — k y ) directions (or vice versa). While the mean field 
study for the case of F 2 (q) = F 2 (q) showed no prefer- 
ence of diagonal order, 11,12 the various experiments in 
cupratcs indicates possibility of both parallel and diago- 
nal fluctuating stripes. 21-23 Here we consider the follow- 
ing model Hamiltonian which offers us to study the inter- 
play between the parallel and diagonal nematic phases. 

ker 

£ [F 2 (q)Ci(k)Ci(k')+^(q)C 2 (k)C2(k')] 

kk'qcrcr' 

xc k+a <J c L-| CT ' c k'+a (T 'C k _q (T , (2) 
where -F 2 (q) and -F 2 (q) are given as follows. 

Here e(k), Ci(k), and C 2 (k) are given by 

e(k) = — t(cos k x + cos k y ) — 2i'cos k x cos k y — /i(4) 
Ci(k) = cos k x — cos k y (5) 
C 2 (k) = 2 sin k x sin k y . (6) 

The mean field Hamiltonian for the uniform nematic or- 
ders is written as 

IAI 2 IA'1 2 

# mean = £ EkC*^ + + h^T , (?) 

where 

e k = e k -ACi(k)-A'C 2 (k). (8) 



Here A and A' measure the strength of the broken x vs. 
y and (x + y) vs. (x — y) symmetries, respectively. 

A = F 2 ^ (cos k x — cos k y ) 6{—e) 

k 

A' = F^^{2sm.k x sm.k y )6{-e) (9) 

k 

We compute the free energy using mean field theory, and 
discuss the phase transition in the following section. 

III. PHASE TRANSITION OF NEMATIC 
ORDERS 

A. Free energy and Nematic order parameters 

The mean field free energy at zero temperature is given 

by 




where 9(e) is step function. Using adaptive 2-dimensional 
integration, 24 we obtain the free energy in terms of chem- 
ical potential and the nematic orders, A and A' Here, we 
ignored the next-nearest hopping t', but we will consider 
it later in sec. Ill D. 

To understand the nature of the transition between 
the diagonal nematic and isotropic phases, let us first set 
F 2 = 0. Fig. 1 shows the free energy as a function of 
the diagonal nematic order, F(A') for several values of 
chemical potential. It is clear that the transition from 
the diagonal nematic phase to the isotropic phase is sec- 
ond order; the diagonal order parameter as a function of 
chemical potential changes continuously. We found that 
A' has finite value only when F 2 exceeds some critical 
value F 2c which depends on the value of A. For the case 
of F 2 — 0, F 2c N is 0.1876 from our numerical calcula- 
tion. Fig. 2 shows that there is no minimum in the free 
energy except A' = for FJ,N < FJ, C N = 0.1876. 

We now turn on F 2 to understand the interplay be- 
tween parallel and diagonal nematic orders. The behav- 
iors of nematic order parameters as one varies F 2 /F 2 are 
obtained by solving the self-consistent equation, (9). All 
possible types of phase diagram for A and A' are summa- 
rized in Fig. 3. As we increase F 2 , the parallel nematic 
order suddenly develops near /i — where the van Hove 
singularity exists. This is illustrated in Fig. 3 (b). The 
suppression of diagonal nematic order due to the devel- 
opment of the parallel nematic order is clearly observed 
as well. However, the region of finite A' does not change 
as long as F' 2 is fixed. As shown in Fig. 3 (c), a fur- 
ther increase of F 2 leads to a wider region of the parallel 
nematic order, which now totally suppresses the diago- 
nal nematic order inside its territory. However, it does 
not take over the whole region of the diagonal nematic 
order. A finite A' region outside the territory of the par- 
allel nematic order is still found. A further increase of 
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FIG. 1: Free energy as a function of A' for various values 
of chemical potential and F^TVo = 0.2077. The value of free 
energy has been rescaled to at A' = 0. 



.002 



/ — - 

o 

W 

/ — - 
■< 
W 




. 1 



A72t 



0.2 



. 3 



FIG. 2: Free energy as a function of A' for several values of 
the interaction, F^ and fi = 0. The free energy has minimum 
at a finite A' only when F2 is larger than a critical value. Ihe 
value of free energy has been rescaled to at A' = 0. 



F2 eventually removes the diagonal nematic phase in the 
picture as shown in Fig. 3 (d). 

It is worthwhile to emphasize the following important 
features of the nematic orders in the coexistence regime. 
The region where the diagonal nematic order is finite is 
always wider than that of the parallel nematic order, if it 
ever exists. As shown in Fig. 4, A has a finite value for 
\fi/2t\ < 0.115, and A' for \fx/2t\ < 0.19. In other words, 
the critical chemical potential for A', /i C 2 is always bigger 
than that for A, fi c i. At /1 = fi c i, A drops to zero discon- 
tinuously and the electron density also changes abruptly, 
which shows the first order phase transition. When A 
goes to zero at /i c i, A' gets enhanced from 0.0679 to 
0.1529 (in unit of 2t). We found that the region of A is 
also shrunk due to finite A', which is further discussed 
in the following subsection. As the chemical potential 
approaches fj, C 2, A' gradually goes to zero reflecting the 



second order phase transition. At this transition, the 
electron density varies continuously and shows only the 
tiny change of its slope. 

We show the behavior of free energies around the criti- 
cal chemical potentials, to highlight these two distinctive 
phase transitions at /i c i and \ici in Fig. 5. Fig. 5 (a) 
shows that the minimum of the free energy for /i/2i = 
occurs at finite A and A'. As the chemical potential 
approaches to the first critical point, |/i c i/2i| = 0.115 
another minimum point begins to develop at A/2t = 
and A'/2t = 0.1529 as shown in Fig. 5 (b). At the critical 
point fi c i, there exist two clear minima as shown in Fig. 5 
(c). The global minimum of (||,^) is changed from 
(0.1035,0.0679) to (0,0.1529) as \fi/2t\ crosses the criti- 
cal point, 0.115. A further increase of fj, finds a unique 
minimum at D = D' = as shown in Fig. 5 (d). 



B. Phenomenological Analysis of two competing 
orders 



To get an insight on competing two nematic orders, we 
analyze the following Ginzburg-Landau (GL) free energy. 



F GL (A,A') 



-A 2 



2 

-—A 2 A' 2 - 



^ + 1 
6 

—A' 4 - 

4 

k 
2 



—A' 6 



(11) 



6 -A 4 A' 2 



-A 2 A' 4 . 



We expand the free energy in terms of the order param- 
eters up to the 6-th orders, since the parallel nematic 
order denoted by A shows the first order transition. We 
set 7' =£3 = for a simplicity, because different values 
of g' and £3 do not affect our qualitative analysis. We 
introduce positive mutual interaction coefficients (£1 > 
and £2 > 0) between two orders, since the two nematic 
orders suppress each other. Here a'(/3) changes its sign 
as it crosses the critical chemical potential /x C 2(/x c i), and 
it is an even function of the chemical potential due to a 
particle-hole symmetry. We consider the following form 
of (3([i) and a'([i). 



0(ji) = /3 (m 2 -a4)- 
a'{/i) = «o(/x 2 - n 2 c2 ). 



1607 



and 



(12) 
(13) 



The order parameters are determined by solving the 
following equations, and shown in Fig. 6 for several values 
of £1 and £2 ■ 



a+ ^A"+ {/3(h) + 2£ 2 A"J A 2 + 7A 4 A (14) 

2" 



= 

= \a'(fi) + 6 A 2 + £ 2 A 4 + /3'A' Z A' 



(15) 



The amplitude of A is not much affected due to the mu- 
tual interaction term. However, the critical chemical po- 
tential, fi c i is shifted in such a way that the region of the 
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FIG. 3: All types of phase diagrams for A and A' as a function of chemical potential obtained by tuning the ratio between two 
interactions F2/-F2, and t' = 0. The solid and dashed lines denote A and A', respectively. See the main text for the discussion 
on (a) -(d). 
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FIG. 4: Nematic orders and electron density vs. chemical 
potential for F 2 N = 0.1 and F 2 'jV = 0.196. The solid and 
dashed lines denote A and A', respectively. The dot-dashed 
line denotes the electron density. The arrows labeled (a) to 
(d) are the values of chemical potentials used in Fig. 5 and 
the discussion in the main text. 
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parallel nematic phase becomes narrower. The modified 
critical chemical potential, fl c i is given by the following 
equation. 



Mel 



I6CV7 



1 + _ 1 ) _ 2 ^ A " 



a 



A) 



(16) 

On the other hand, the magnitude of A' is suppressed 
due to finite A, and it is modified as follows. 



A' = 



c*o(/4-|» 2 )-£iA 2 -6A 4 



0' 



(17) 



The critical chemical potential, /x C 2 is not affected as long 
as /i C 2 > p-d- Our GL analysis captures the key features 
of our results presented in sec. Ill A. They suppress each 
other in qualitatively different ways. The parallel ne- 
matic order suppresses the amplitude of the diagonal ne- 
matic order, while the window of the diagonal nematic 
order is not affected when they coexist. On the other 
hand, the diagonal nematic order shrinks the parallel ne- 
matic order region, but hardly suppresses the amplitude 
of A. ' 



FIG. 5: The density plot of free energy as functions of A and 
A' for several values of chemical potentials indicated by the 
arrows in Fig. 4, and for given values of interactions, F2N0 = 
0.1 and F2N0 = 0.196. Darker is lower free energy. The 
minimum points for each figure denote the equilibrium values 
of A and A'. The discussion on (a) - (d) can be found in the 
main text. 



C. Fermi surface and density of states 

To understand the nature of the phase transition, we 
investigate effects of nematic orders on the density of 
states (DOS) and Fermi surface. The DOS for several 
values of A and A' is shown in in Fig. 7. Without the 
nematic orders, DOS has a singularity at E + \i = 
originated from the van Hove singularity (VHS). It was 
shown that the development of A (the dotted line) leads 
to a dramatic change in DOS. The parallel nematic or- 
der splits the VHS into two peaks occurring near the van 
Hove filling. As a result, the free energy is lowered. 12 
This feature is inherited from the logarithmic singular- 
ity in the free energy. (See ref.12 for detail.) On the 
other hand, A' is nothing to do with the VHS as shown 
in Fig. 7. The peak from VHS is not affected by a de- 
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FIG. 6: Order parameters determined from the Ginzburg- 
Landau equation of (11) for various values of £i and £2. The 
other parameters are a = 0.1, j3o = 4.0, 7 = 1.0, a' Q = 1.0, 
P' = 1.0, 7' = 0.0, £3 = 0.0, fid = 0.25, and fj, c2 = 0.5. The 
solid and dashed lines denote A and A', respectively. 
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FIG. 7: Density of states for various values of A and A'. 
The solid line is for (A/(2t), A'/(2t)) = (0,0), the dashed 
for (0,0.25), the dotted for (0.25,0), and the dot-dashed for 
(0.25,0.25). A finite A splits the van Hove singularity at 
E = — into two peaks, while A' shows the minor effect on 
the density of states. 



vclopment of A'. The DOS has only minor change due 
to A', which is a slight enhancement near the band edge 
and suppression near the center of the band. Therefore, 
the free energy develops a minimum continuously from 
A' = to a finite A', as one changes fi, thus the transi- 
tion between isotropic and the diagonal nematic phases 
is second order. 

One of important consequences of the nematic order is 
a deformation of the Fermi surface. The deformation of 
Fermi surfaces for various values of A and A' are shown 
in Fig. 8. Fig. 8 (a) shows the undeformed Fermi sur- 
face to make comparison with (b)-(d). A finite paral- 
lel nematic order, A squeezes the Fermi surface along a 
parallel axis of the lattice as shown in Fig. 8 (b). For 
example, it shrinks the Fermi surface in fc^-direction and 





(c) 4=0.0, A'=0.4 



(d) A = 0,4, A 





FIG. 8: Fermi surface for several values of A and A'. The 
deformation of the original Fermi surface shown in (a) along 
the parallel (diagonal) axes of the lattice is due to the parallel 
(diagonal) nematic order. Examples of Fermi surface for the 
parallel and diagonal nematic order phases shown in (b) and 
(c), respectively. An example of Fermi surface for finite A 
and A' is shown in (d). 



expands it in fc y -direction, or reverse way. On the other 
hand, the diagonal nematic order A' deforms the Fermi 
surface along a diagonal axis of the lattice as shown in 
Fig. 8 (c). For example, it shrinks the Fermi surface in 
(k x — /cj,)-direction and extend in (k x + /cj,)-direction, and 
vise versa. The discontinuous change of A related to the 
van Hove singularity leads to a dramatic change in the 
shape of Fermi surface as shown in Fig. 8 (b). On the 
other hand, the deformation along the diagonal direction 
develops continuously, as one changes /z. It is important 
to note that D' does not affect four points on the Fermi 
surface, (±kp Xl 0) and (0, ±kp y ), which eventually lead 
to the van Hove singularity at the van Hove filling and 
the formation of the parallel nematic phase inside the 
diagonal nematic phase. 



D. Effects of next- nearest hopping 

In this section, we introduce the next-nearest hopping 
integral, t' in (4), which breaks a particle-hole symmetry. 
We found that a negative sign of t' shifts the region of ne- 
matic phases toward hole-doped region. However, since 
the parallel nematic order always appears near a VHS, 
while the diagonal nematic order is not directly affected 
by a VHS, we expect that the region of parallel nematic 
and diagonal nematic phases will be shifted in a slightly 
different way, as we increase t' . On the other hand, qual- 
itative features such as nature of phase transition are not 
affected by a finite t' . 
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FIG. 9: Nematic orders and electron density vs. chemical 
potential for t' = -OAt, F 2 N = 0.1 and F 2 'iVo = 0.196. The 
solid and dashed lines denote A and A', respectively. The 
dot-dashed line indicates the electron density. 



Fig. 9 shows a phase diagram for t' — — OAt. The par- 
allel nematic order is shifted more to the hole doped re- 
gion than the diagonal nematic order. At the half-filling, 
only the diagonal nematic order is finite. As we increase 
(or hole concentration), the coexistence of two ne- 
matic phases appears. A further increase of leads to 
a suppression of the diagonal nematic order, while the 
parallel nematic phase gets stronger. Finally only the 
parallel nematic order remains, which eventually disap- 
pears in the phase diagram, as the hole concentration is 
increased. The nematic order rotates from the parallel 
order to the diagonal order, as we increase hole doping 
concentration, and two phases coexist near the bound- 
ary. This feature resembles the rotation of charge fluc- 
tuations associated with Bi g channel to B\ g channel by 
increasing doping concentration which was reported in 
the Raman spectroscopy measurement on I^-zSr^CuCU 
(LSCO). 23 The rotation of spin modulation by changing 
doping concentration in LSCO was also found in elastic 
neutron scattering patterns. 22 Elastic neutron scattering 
studies on LSCO showed one-dimensional spin modula- 
tion along the orthorhomic 6-axis in lower doping con- 
centrations, and another type of spin modulation paral- 
lel to the tetragonal axes in high doping concentrations. 
The coexistence of two types of spin modulations near 
the boundary was also reported. 22 While direct compar- 
isons to the neutron scattering patterns and/or Raman 
spectroscopy data require further theoretical studies on 
corresponding susceptibilities in the nematic phases, we 
expect that the behavior of Fermi surface deformation 
within our model offers a consistent picture compared 
with the experimental observations. 

Fig. 10 summarizes the typical types of phase diagrams 
for two nematic orders with the next-nearest hopping. 
For a finite negative t' , both A(/i) and A'(/i) are shifted 
to hole-doped region but in a slightly different way, as 
we discussed. When they coexist, the maxima of two 
nematic phases do not coincide as shown in Fig. 3 due to 



their different t' dependences. However, the qualitative 
behaviors of competition between two nematic orders are 
similar to the case without t' presented in the sec. Ill A. 



IV. DISCUSSION AND SUMMARY 

A quantum analog of classical liquid crystal in terms 
of broken symmetry has been discussed in the context 
of a doped Mott insulator. 2 Among the electronic liquid 
crystal phases, nematic phases can be viewed as fluctu- 
ating stripes whose segments are fluctuating in time but 
oriented in a particular direction, and hence breaks ori- 
cntational symmetry. 

Recent neutron scattering measurements of detwinned 
YBa2Cu307_,j have indicated a possible existence of two 
dimensional anisotropic liquid crystalline phase in high 
temperature cuprates. 25-27 Extensive neutron scattering 
measurements of La2-xSr x Cu04 in a wide range of dop- 
ing have also revealed the doping dependence of the static 
or quasi-static spin ordering in insulating and supercon- 
ducting phase. 21,22 An interesting observation is that the 
orientation of the spin modulation depends on the dop- 
ing concentration. It was found that the spin modula- 
tion vector is diagonal to the Cu-0 bond in the insu- 
lating spin glass phase, while inside the superconducting 
phase it is parallel to the tetragonal axes. These two 
types of spin modulation coexist near the boundary be- 
tween the insulating and superconducting phases. Such a 
one-dimensional nature of the spin correlations is consis- 
tent with a stripe-like ordering of the holes in the Cu02 
planes. 3 

Inelastic light-scattering spectra of underdoped 
La2- x Sr x Cu04 single crystals showed additional Drude- 
like responses in Bi g and B\ g channels for x — 0.02 
to x = 0.10, respectively. 23 This was interpreted as 
experimental evidence of fluctuating charge stripes 
whose orientation rotates from diagonal to parallel by 
varying the doping concentration of Sr from x = 0.02 
to x = 0.10. 23 More extensive studies for various 
doping concentrations will be required to determine the 
coexistence of two types of charge modulations. While 
direct comparisons to these experimental data require 
theoretical studies on corresponding susceptibilities, 
we expect that our phenomenological model provides 
a possible explanation of the rotation of the spin and 
charge modulations observed in LSCO. 

In summary, we have studied the interplay between 
parallel and diagonal nematic phases using the phe- 
nomenological model Hamiltonian within mean field the- 
ory. We found that the parallel and diagonal nematic 
phases compete each other - the parallel nematic order 
suppresses the amplitude of the diagonal nematic order, 
while the diagonal nematic order shrinks the window of 
the parallel nematic order. However, they still coexist in 
a finite window of parameter space. The transition to the 
parallel nematic phase is strongly first order, while the 
diagonal nematic phases shows the continuous transition 
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FIG. 10: Various types of phase diagrams for A and A' with the next-nearest-hopping t'. The solid line is designated for 
A and the dashed line for A'. A negative t' shifts both A and A' toward hole-doped region, (a) F-i — and F'^ > 0. (b) 
A small increase of F2 leads to a development of the parallel nematic phase inside the diagonal nematic phase region, (c) A 
further increase of F2 increases the parallel nematic phase region, thus suppresses the region of diagonal nematic phase, (d) 
the diagonal nematic order is totally suppressed within the region of parallel nematic phase, while it still survives outside the 
region, (e) The parallel nematic order eventually removes the diagonal nematic order by expanding its region. 



to the isotropic phase. Effects of quantum fluctuation of 
the diagonal nematic order parameter on various quan- 
tities near a quantum critical point are the subjects of 
future studies, and the single particle self-energy correc- 
tion due to the fluctuation of a collective mode will be 
presented elsewhere. 28 
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